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We investigate the metal insulator transitions at finite temperature for the Hubbard model with 
diagonal alloy disorder. We solve the dynamical mean field theory equations with the non crossing 
approximation and we use the coherent potential approximation to handle disorder. The excitation 
spectrum is given for various correlation strength U and disorder. Two successive metal insulator 
transitions are observed at integer filling values as U is increased. An important selective transfer 
of spectral weight arises upon doping. The strong influence of the temperature on the low energy 
dynamics is studied in details. 



I. INTRODUCTION 

Strongly correlated materials have been the subject of 
intensive researches for the last forty years. The simplest 
lattice model accounting for correlations is the one band 
Hubbard model. Great progresses have been achieved 
within the last ten years with the development of the 
dynamical mean field theory (DMFT) (for a review see 
Ref. Q). One important result of the DMFT is the de- 
scription of the metal to insulator transition driven by 
correlations (Mott transition 2 ) for the half- filled Hubbard 
model. Today, the Mott metal insulator transition (MIT) 
is the subject of intensive experimental and theoretical 
studies^ Recently, a renew of interest for the Mott tran- 
sition has occurred in relation with the so-called orbital 
selective Mott transitions*^ observed in multi-orbital 
compounds. 

In most of the experimental studies, strong electron 
correlations are closely connected with disorder effects. 
So, to change experimentally the band-filling, the band- 
width or other parameters, one replaces usually one 
atom of a given compound by another one. That is 
true for strongly correlated transition metal oxides or 
similar compounds. For example, replacing La by Sr 
in La^Sri-^TiOs changes the doping, i.e. the band- 
filling&a In Ca^^Sr^ V0 3f i2iii Ca and Sr atoms have the 
same valence but different atomic radii, which changes 
the angle of the V-O-V bond and, correspondingly, the 
band-width. More recently, a very interesting studyi^ on 
the metal-insulator transition in the disordered and cor- 
related system SrTii-^Ru^Oa has been performed. From 
their transport and optical data, the authors propose a 
classification scheme with six kinds of electronic states 
depending on x. The disorder and the electron correla- 
tion effects should be considered together to understand 
the measured electronic structure evolutions. 

A generic model to study the common influence of dis- 
order and correlations is the Hubbard model including di- 
agonal disorder. As we will outline below, it shows many 



features and metal insulator transitions not present in the 
standard Hubbard model. A possible method to solve the 
disordered Hubbard model is the combination of DMFT 
with the coherent potential approximation (CPA)i4 to 
treat disorder. It was used in the seminal work of Laad 
et al.r^ where the local impurity model of the DMFT 
was solved by the iteration perturbation theory (IPT). 
Interesting results concerning the stability of the Fermi 
liquid metal against small disorder and the effect of the 
interplay between correlations and disorder on the optical 
and Raman response had been obtained. 

For the disordered Hubbard model, a first classification 
of the metal-insulator transitions at non-integer fillings 
(and T = 0) was given in Ref. [l5|. It used an analogy to 
the Zaanen-Sawatzky- Allen scheme of charge-transfer or 
Mott-Hubbard insulators, respectively (see also Ref. 
for the relationship to Anderson localization). Here, we 
will concentrate first on the Hubbard model with diag- 
onal alloy disorder at integer filling. With increasing 
Hubbard correlation U we find two metal-insulator tran- 
sitions: first from a non-correlated band insulator to a 
metal and a second transition from a metal to a Mott- 
Hubbard insulator. 

The transition from a band insulator to a Mott- 
Hubbard insulator has been the subject of intensive cur- 
rent studies within the framework the half-filled ionic 
Hubbard model (IHM).iL±& The IHM corresponds to the 
spatially ordered limit of the binary alloy system which 
is considered in the present paper, A X B\- X , where the 
A atoms are separated from each-other by an identical 
sequence of B atoms. In particular, intensive analyt- 
ical and numerical efforts have been performed to an- 
alyze the band insulator (BI) to Mott insulator (MI) 
transition in the Hubbard chain with alternating on-site 
energiesi 19 i 20 i 21 i 22 i 23 i 24 These studies show that the half- 
filled ionic Hubbard chain has two transitions as U is 
increased. The first is from the BI to a bond-ordered, 
spontaneously dimerized, ferroelectric insulator (FI). The 
second, when U is further increased, involves a transition 
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between the FI and the MI. In this phase diagram the 
very unconventional metallic state is realized only at the 
BI-FI transition point. 

Very recent studies of the generalized AB n -i IH chain 
for n > 3, where n — 1 is the length of domain of _B-atoms 
which separates A atoms, show that at commensurate 
band-filling l/n the system shows, with increasing t/, 
a similar sequence of transitions BI-FI-MI>2^ Therefore, 
in the case of ID ordered binary chains, the BI phase is 
always separated from a MI by the insulating ferroelectric 
phase. Although the Mott transition in the disordered ID 
binary chain has not been studied in detail, the Anderson 
localization 26 excludes the presence of a metallic phase 
in a finite sector of the ground state phase diagram of a 
ID disordered IH model. 

Therefore, the approach presented in this paper al- 
lows us to investigate in detail the band-insulator to 
correlated metal and metal to Mott-Hubbard insulator 
transitions in a disordered binary-alloy system for higher 
dimensions,? 7 Below we analyze these metal-insulator 
transitions in detail. For that purpose we consider the 
spectral weight transfer with doping, and we find that it 
occurs in a selective way, different for the two disorder 
sites in question. 

In the present work, we propose a finite temperature 
calculation of the dynamical properties by using the non- 
crossing approximation (NCA)SSi2S to solve the DMFT 
local Anderson problem in the self-consistent loop. Then 
we combine the DMFT with CPA which is a very natural 
way since both methods use a local self-energy. Spectral 
densities and spectral weight transfers will be computed 
with the NCA. The role played by the temperature for 
the low energy physics will be discussed in details. 

In the following section the model Hamiltonian is pre- 
sented. This section also presents the application of the 
DMFT on this Hamiltonian and a solving method based 
on the NCA. In the third section, the main results are 
discussed. These are the insulator-metal-insulator tran- 
sition at integer filling and a detailed analysis of possible 
transitions at integer and non-integer fillings by using 
their spectral properties. 



II. HAMILTONIAN AND SOLVING METHOD 

We consider the following Hamiltonian, where correla- 
tions and diagonal disorder are present : 

H= -t + H.c) 

<i ,j>,iy 



(1) 



In expression the sum < i, j > is the sum over nearest 
neighbor sites of a Bethe lattice. cf a (respectively do) de- 
notes the creation (respectively annihilation) operator of 
an electron at the lattice site i with spin a and Ui a is the 
occupation number per spin. U is the on-site Coulomb 



repulsion and t is the hopping term. A fraction 1 — x of 
sites (sites A) have a local on-site energy £4 — ea and 
a fraction x of sites (sites B) have a local on-site energy 
£j = £b- We will note A = ea — sb the energy difference 
between the two types of sites. 

Following the DMFT procedure of Ref. [H by inte- 
grating out all fermionic degrees of freedom except for a 
central site i — o, the lattice model Q can be mapped 
onto an effective impurity model. The expression of the 

A B 

corresponding Hamiltonian H ct j depends on the nature 

A B 

of the central site i = o. It contains a local part H, ' 
and a part corresponding to the coupling to the effective 
medium. This coupling is determined self-consistcntly. 
We then write : 



where 



t a,b _ u a,b tt 



H tf = Un Q ^n ol 



(2) 



is depending on the nature (A or B) of the central site i = 

A B 

o. In the expression of -Hi' , we have written explicitly 
the chemical potential fi. The coupling Hamiltonian with 
the effective medium is 



ffmcd = ( W k b t C °° + HX 

ka 

+ V Etbi b r . 



Wt represents the hybridization between the site i = o 
and the effective medium, is the band energy of the 
effective medium, bi (respectively b^) is the creation 
(respectively annihilation) operator of an electron in the 
effective medium. The effective medium can be charac- 
terized by the effective dynamical hybridizations : 



J(u>) 



l^l 2 



^ uj + i0+ - £,- 
k 



On a Bethe lattice, the self-consistent equations of the 
DMFT can be simply written 

J{Lo)=t 2 G a {uj) , 

where G a (u>) is obtained by averaging over disorder with 
the CPA procedure. We have G a (w) = xG^(lu) + 
(1 - x)Gf (u>) where G$(u) (respectively G^{u>)) is the 
Green's function solution of the effective local impurity 
problem represented by the Hamiltonian (respec- 
tively 

To solve the local impurity problem H^ c (where a = 
A, B) with the NCA, we introduce four local states \m a ) 
and their local energies E^. m = 1 corresponds to an 
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empty local site, m = 2 (respectively m = 3) corre- 
sponds to a singly occupied local site with an electron 
of spin | (respectively j), an( l m — 4 corresponds to a 
doubly occupied local site. Of course, energy E" is a 
independent and because of spin degeneracy, we have 
E% = Eg = e a — /i. For the doubly occupied site 
E% = 2e a - 2/i + U. 

The NCA equations^ read : 

SfM = 2 J dep(e)f(e)P 2 a (u + e) 

= y dep(e)f{-e)P? (w - e) 

+ /" <fep(e)/( £ )P 4 Q (w + £ ) 
E£(w) = ^(w) 

Ej(w) - 2 /" rfe ( 5( £ )/(-e)P 2 Q (a ) -e) 
where the local propagators P^ are defined by : 



^ 1 c + zO+ 

and the effective medium spectral density by p(e) = 
— -Im(J^(uj)). f is the Fermi-Dirac distribution at tem- 
perature T. We finally obtain the site specific one particle 
Green's function : 

G«H = -L Jdee-f > °VS(e)P?(u> + e)+fi(e)P?(w+e) 
- P Ue)P Q a (e ~ «)* - P a 2 {e)P?{e - «)*] 

where pjjjfe) = — -Im(P^(e)) and the partition function 
is Z$ = Jdee-f >e frS(E) + 2pf(e) +p$(e)] . 

III. RESULTS AND DISCUSSION 

In this section, we present the calculated densities of 
states for the disordered Hubbard model. In the first sub- 
section, we focus on the possible metal insulator transi- 
tions and in particular the nature of these transitions. 

In the second sub-section MI Bl we show how the site 
selective spectral weight transfer can help to understand 
the interplay of disorder and correlations for dynamical 
properties. 

A. Metal insulator transitions 

One has to distinguish the metal insulator transitions 
at integer filling (n = 1) from those which occur at non- 
integer fillings (n = x or n = 1 — x) since at A = the 
MI transition with increasing on-site repulsion U takes 
place only in the former case. 

Let us first concentrate on those with n = 1. In that 
case, the concentration of sites A has to be x — 0.5. 



Then, for each A and sufficiently large values of U we 
observe an insulating state for n — 1. Before discussing 
the commensurate n = 1 case it is very instructive to 
consider the case of a doping slightly below half-filling. 
In this case we have a strongly correlated metal with 
heavy quasi particles for all U. That is illustrated in Fig. 
□ for A = 2 eV, U = 6 eV, and n = 0.9. 

NCA applied to DMFT is known 29 to supply a fairly 
good solution even away from half-filling. For this 
strongly correlated situation, we found (plain line) a 
metallic state with four incoherent bands and one co- 
herent resonance at the Fermi level, witch is a signature 
of the presence of quasiparticles. The temperature de- 
pendence of the quasiparticle peak is shown in part (b) 
of the figure. Although the density of states is almost 
constant at the Fermi level, the low energy states near 
the Fermi level are strongly temperature dependent. In 
part (a) of figure we plot as well the spectral densities 
obtained from the solutions of the two impurity mod- 
els with a site A or B embedded in the same effective 
medium. This gives enlightenment on the nature of the 
different bands forming the density of states for the dis- 
ordered system and then allows to clarify the nature of 
the near insulating state. Using the NCA, it is possible to 
reduce arbitrarily the doping with increasing numerical 
efforts. Band Al (respectively Bl) is mainly composed of 
singly occupied A sites (respectively B sites) and spectral 
weights are proportional to the probability of removing 
one electron from a site A (respectively B). Correspond- 
ingly, the spectral weight of band A2 (respectively B2) 
is proportional to the probability of adding one electron 
to a singly occupied site A (respectively B). The energy 
gap between band Al and band A2 (or Bl and B2) is 
of the order of U = 6 eV. We can identify Al to a so 
called lower Hubbard band (LHB) and B2 to a so called 
upper Hubbard band (UHB). We can therefore consider 
that the system is a strongly correlated metal, close to a 
Mott insulator, with an insulating gap between the LHB 
(Al) and the UHB (B2) and a reduced effective value 
of the correlation strength U — A. Note the asymmetric 
role played by the two types of site. For n — 0.9, charge 
fluctuations are already blocked by correlations on sites 
B which are half-filled. The Fermi level resonance peak 
is therefore only built up by electrons evolving on sites 
A. This result is confirmed by the respective local occu- 
pations ha — 0.4 and tib = 0.5. 

Part (a) of Fig. shows the densities of states for a 
decreasing U (U = 7, 5, 4, 3, 2 eV) and a given A = 4 eV 
at T = 1000 K for x = 0.5 and n = 1. We observe 
two phase transitions. As expected, starting from the 
strongly correlated situation U = 7 eV, the system un- 
dergoes a transition to a metallic state for decreasing U . 
There is a second critical value at which the system en- 
ters into a new insulating state with further decreasing U . 
The nature of the insulator is completely different here. 
This can be clarified by examining part (b) of figure ^ 
Sites B are almost filled (ub = 0.915) whereas sites A are 
almost empty {ua — 0.045). The two side bands of the 
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FIG. 1: (a) Density of states (plain line) obtained from 
DMFT-NCA and CPA for U = 6 eV, A = 2 eV and 
T = 1000 K. Dotted line (respectively dashed line) shows the 
corresponding density of states for the impurity problem with 
a site A (respectively B) surrounded by the effective medium, 
(b) Temperature dependence of the low energy part of the 
calculated spectrum. 



gap correspond to singly occupied sites and cannot be 
interpreted in terms of Hubbard bands. Therefore, this 
insulator can be understood as an ordinary band insula- 
tor. The metallic state for intermediate values of U can 
be understood since in that case the energies for doubly 
occupied B sites eb + U and singly occupied A sites ea 
are nearly degenerate. For U = A we have a special sit- 
uation (see Fig. . In that case we have a single peak 
in the DOS at the Fermi energy, whereas in all other 
situations a double peak develops. 

We believe that the series of phase transitions 
insulator- metal-insulator is a rather generic situation and 
should occur in a wide class of models with two different 
species at equal concentration in dimensions higher then 
one. Very recently a similar insulator-metal-insulator 
transition was obtained for the ordered Hubbard model 
with two sites on a bipartite latticed Also, there is a 
quite large similarity to the transition between a Mott- 
Hubbard insulator and a band-insulator observed in the 
ID ionic Hubbard model, however due to the peculiarity 
of ID systems the metallic state in this case is reduced 
to only one critical point corresponding to the transition 
from a BI to FI. 

To characterize the different phases more in detail we 
computed also the effective masses in a slightly doped 
situation (n — 0.96, not shown). We used the equation 
m*/m = 1 — [3Re (E"(cj)) /du] u=0 . Crossing the upper 
insulator-metal transition, we found a strong reduction 
of m* B /m from 13.5 (for U = 7 eV) to 2.4 (for U = 5 eV) 
while m* A /m evolves slightly from 1.5 to 1.9. This is 
in contrast to the band-insulator case (U = 2 eV) with 
the effective masses m* A /m — 1.004 and m B /m = 1.07 




FIG. 2: (a) DMFT-NCA-CPA densities of states for U = 
7, 5, 4, 3, 2 eV, A = 4 eV, T = 1000 K, x = 0.5 and n = 1.0. 
(b) Corresponding results for site selective densities of states. 



proving the absence of heavy quasiparticles. 

Finally, let us discuss the metal-insulator transitions at 
non-integer fillings. They have been found by Byczuk et 
alm^ Using the numerical renormalization group method 
at zero temperature, they have shown that the system 
becomes a Mott insulator at strong interaction for n = x 
or n — 1 — x. Here we propose a visual interpretation of 
this result. Figure |3| displays the densities of states for 
U = 2A = 4 eV, n = 0.9 and various disorder concentra- 
tions x of site B between and 1. For x = 0, there are 
only sites A whose spectral densities are distributed over 
two Hubbard bands, and the lower one is partially filled. 
Increasing the concentration of B sites leads to an in- 
creasing strength of the two Hubbard peaks correspond- 
ing to B (Bl and B2). For x = 0.9, the lower Hubbard 
peak of site B is completely filled and sites A are nearly 
empty, leading to an insulating situation. To observe the 
insulating state at x = 0.1 we should exchange the role 
of A and B site (A < 0). 



B. Selective spectral weight transfer 

In this section, we investigate a larger range of doping 
values, from the quarter-filled situation where n = 0.5 to 
n = 1.5, but keeping x — 0.5. We found an important 
transfer of spectral weight and explain it by using the se- 
lective spectral densities supplied by our approach. Two 
particular cases of non-integer filling n = 0.5 and n = 1.5 
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FIG. 3: Density of states obtained from DMFT-NCA and 
CPA for U = 2A = 4 eV, n = 0.9 and various disorder 
concentrations x of site B between and 1. 








' \ Al 


B2 If " 


/ ^2 


BlV 



energy [eV] 




5 
energy [eV] 



FIG. 4: (a) Densities of states for U = A = 4 eV, T = 1000 K 
and for various fillings from n = 0.55 ton = 1.45. (b) Spectral 
weights of the different bands, (c) Corresponding densities of 
states for sites A. (d) Corresponding densities of states for 
sites B. The position of the Fermi level is marked with an 
asterisk on each spectrum. 



IV. CONCLUSION 



will be discussed in details, with a new type of insulating 
state. In part (a) of figure 0] the densities of states are 
displayed for U = A = 4 eV, T = 1000 K and for vari- 
ous fillings from n = 0.55 to n = 1.45. The position of 
the Fermi level is marked with an asterisk on each spec- 
trum. As the filling is increased, an important transfer of 
spectral weight takes place between the different bands. 
The spectral weights are reported in part (b) of the fig- 
ure. Part (c) (respectively part (d)) of the figure shows 
site selective densities of states for sites A (respectively 
B). We found that the observed spectral weight trans- 
fer is characteristic of strongly correlated bands and is 
site specific. For n going from 0.55 to 0.95 the spectral 
weight transfer occurs for sites A (from 1.05 to 1.45 for 
sites B). For n = 0.55 and n = 1.45 the system is close 
to an insulating state. Nevertheless, the nature of the 
insulating phase is strongly different from the half-filled 
situation, where we have a Mott insulator and two half- 
filled spectral densities for sites A and B. For the almost 
quarter-filled situation (n = 0.55) we found an hybrid 
type of insulating phase, where sites B are half-filled and 
sites A are almost empty. For B sites, charge fluctua- 
tions are freezed by correlations like in the usual Mott 
insulator. For A sites, the blocked microscopic processes 
are of charge transfer type like in a band insulator. For 
n = 1.45 we have the symmetric situation. 



In this paper, we proposed an approach based on the 
dynamical mean field theory which is able to handle diag- 
onal disorder in a strongly correlated Hubbard model at 
finite temperature for any doping. We used the coherent 
potential approximation in the self-consistent equation of 
the DMFT and we applied the non-crossing approxima- 
tion to the local impurity problem on which the lattice 
model is mapped. 

For a decreasing U we showed that the disordered half- 
filled system undergoes two successive metal insulator 
transitions of different nature. The first transition is a 
Mott like transition with a reduced effective correlation 
strength. The second transition is a charge transfer like 
transition. We pointed out that this situation is a very 
generic one. We also discussed the differences and sim- 
ilarities between the considered disordered binary-alloy 
model and the ionic Hubbard model in one and higher 
spatial dimensions. 

In addition to the insulating phase at half filling where 
both types of site are half-filled and all charge fluctu- 
ations are blocked by correlations effects, we found an 
hybrid type of insulating state away from half filling, i.e. 
at n = x or n — 1 — x. We illustrated it for the quarter 
filled and the three-quarter filled situations. Charge fluc- 
tuations are blocked by correlations for one type of sites 
and by charge-transfer excitations for the other type of 
sites. 
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It could be interesting to include in the present model 
two orbitals with different bandwidth to investigate the 
role of disorder in the orbital selective Mott transition. 
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